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PROPAGATION OF A FREE FLAME IN A TURBULENT GAS STREAM 

By William R. Mickelsen and Norman E. Ernstein 


SUMMARY 

Elective flame speeds of free turbulent flames were measured 
by photographic, ionization-gap, and photomultiplier-tube 
methods, and were found to have a statistical distribution 
attributed to the nature of the turbulent field. The effective 
turbulent flame speeds for the free flame were less than those 
previously measured for flames stabilized on nozzle burners, 
Bunsen burners, and blufl bodies. The statistical spread of 
the effective turbulent flame speeds was markedly wider in the 
lean and rich fuel-air-ratio regions, which might be attributed 
to the greater sensitivity of laminar flame speed to flame temper- 
ature in those regions. Values calculated from the turbulent 
free-flame-speed analysis proposed by Tucker apparently form 
upper limits for the statistical spread of free-flame-speed data. 
Hot-wire anemometer measurements of the longitudinal velocity 
fluctuation intensity and longitudinal correlation coefficient 
were made and were employed in the comparison of data and m 
the theoretical calculation of turbulent flame speed. 

INTRODUCTION 

The high volumetric heat-release rate required in present 
and future jet-engine combustors has placed an ever-growing 
emphasis on turbulent combustion research. The theory of 
turbulent flame propagation is as yet in a formative stage, 
partly because of the lack of reliable experimental methods 
for determining the effect of turbulence on the rate of flame 
propagation. The purpose of this investigation was to ex- 
plore a new experimental method believed to approximate 
more nearly a theoretical, or ideal, flame model than previous 
methods. 

The major portion of experimental data on turbulent flame 
speeds has been taken in open flames stabilized on Bunsen 
burners (refs. 1 to 3) or on Mache type burners (refs. 4 and 5). 
A substantial amount of experimental data has also been 
obtained from flames confined in a duct and stabilized on 
bluff bodies (refs. 6 to 8). Another experimental method 
consists in measuring the speed of a flame advancing into a 
turbulent fuel-air mixture flowing through a tube (ref. 9). 
In general, the definition of turbulent flame speed used in 
these methods is that given in reference 8: 

S T =U i? (!) 

A 

(All symbols arc defined in appendix A.) Equation (1) has 
been used for local measurements by the flame-front angle 
method and for over-all flame-speed measurements b} the 


flame-area method. These two methods are shown in the 
folio whig sketches: 

' A" J 

„ _,.A „ 

S T ~ U sin a 



U t— 

Flame- front angle method 



Flame- area method 


Turbulent-flame-speed values obtained by these methods 
may be significantly affected by conditions external to the 

flame front itself, such as . 

(1) Large fluctuations in instantaneous flame-front, posi- 
tion that introduce uncertainty in the determination of 
mean position in long-time-exposure photographs of 
stabilized flames (ref. 10) 

(2) Existence of piloting zones at the runs of Bunsen 
burners, immediately downstream of bluff bodies, or in 
the boundary layer of flame tubes, which may affect the 
burning rate of the flames near such zones (refs. 7 and 

(3) 1 Large velocity differences between unburned and 
burned gas flows, which may introduce considerable 
turbulence in stabilized flames (refs. 2, 6, 7, and 11) 

(4) Curvature of unburned-gas stream lines at the flame 
front which mav introduce substantial error in flame 
speeds measured in stabilized flames by the angle 
method (refs. 7 and 12) 

(5) Variations in local flame speeds along the flame enve- 
lope, which may introduce substantial error in average 
flame speeds measured by the flame-area method 

This investigation was conducted to determine the effect 
of turbulence on the growth of a free flame. It was thought 
that the free flame would not be subject to the complicating 
conditions present in stabilized flames such as those listed. 
In addition, the free-flame growth could take place entirely 
in a homogeneous, isotropic, turbulent field away from 
boundary layers and bluff-body wakes, thus providing a 
more nearly fundamental flame model. 

The method used in the present investigation is an exten- 
sion of the “soap-bubble” technique (ref. 13) and consists 
of the observation of the growth of a free-flame globule as 
it is carried downstream in a flowing, turbulent, homo- 
geneous mixture of gaseous fuel and air. The flame globule 
was initiated by a single spark, and its growth wa.s recorded 
by three separate methods, high-speed motion-picture 
photography, ionization-gap probes, and photomultiplier- 
tube signals. 


1 Supersedes NACA TN 3456, 


“Propagation of a Free Flame in a Turbulent Gas Stream 


*> by William R. Mickelscfi and Norman E. Ernstein, 1955. 
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Tile flame-globule growth was characterized by an effec- 
tive turbulent flame speed St, which is defined by a mass- 
balance equation similar to that derived in reference 13: 


p F dr 
Pa dt 


( 2 ) 


by the use of three interchangeable grids, and the param- 
eters of the turbulence were measured with constant- 
temperature hot-wire-anemometer instrumentation. 

I he investigation was carried out at the Lewis laboratory 
as a part of the combustion research program. 


Equation (2) can be clarified by the following sketch : 



The rate of change of the globule radius was measured over 
periods of time less than the characteristic time of the tur- 
bulence. This means that, on the avei’age, the flame front 
traveled through less than one turbulent eddy during the 
time of observation. Tn the past, average turbulent flame 
speeds have generally been measured from photographs 
having long exposure times compared with the characteristic 
time of the turbulence. 

The experimental work was carried out in two phases. 
The first phase was initiated and carried out at the NACA 
Lewis laboratory and is reported in reference 14. (Because 
of its limited circulation, considerable portions of the work 
reported in ref. 14 are included in this report.) Turbulent- 
flame-speed measurements were obtained from high-speed 
motion-picture records of the growth of the free-flame glob- 
ule as it was swept downstream in an enclosed tunnel. 
The turbulent field was introduced by means of a wire 
grid placed at the tunnel inlet some distance upstream of 
the spark-electrode position. The turbulent field was 
characterized by hot-wire-anemometer measurements. 

The effective turbulent flame speeds determined by the 
photographic method varied for runs made under identical 
fuel-air and mean stream conditions. The variation was 
concluded to be due to the statistical nature of the turbu- 
lent field. Statistical analysis of groups of 30 or more runs 
made at identical fuel-air and mean stream conditions 
showed that the data for each group approximated a normal 
probability distribution. When plotted on the basis of 
cumulative probability of occurrence, the effective turbulent- 
flame-speed data showed a consistent increase with increasing 
turbulence intensity. 

The second phase of the investigation was carried out in 
a free jet with a technique whereby the cumulative proba- 
bility of occurrence of turbulent flame speed could be deter- 
mined from flame groups consisting of thousands of separate 
flame globules. This technique utilized ionization-gap 
probes and counters, which indicated directly the percent- 
age of flames reaching or exceeding any given diameter at 
a series of stations downstream of the spark electrodes. 
By this method, flame speeds were measured over a range 
of fuel-air weight ratios from 0.053 to 0.090 and a range of 
mean stream velocities from 35 to 142 feet per second at 
a stream static pressure of 1 atmosphere and a stream static 
temperature of 85° F. The stream turbulence was varied 


APPARATUS 

AIR-DUCTING. WIND-TUNNEL, AND AIRFLOW INSTRUMENTATION 

The air-ducting, wind-tunnel, and airflow instrumenta- 
tion used in the first phase of the investigation and reported 
in reference 14 is shown in figure 1(a). Alterations made to 
the apparatus for the second phase of the investigation are 
shown in figure 1(b). The alterations essentially converted 
the system from an enclosed tunnel to a free jet. 
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Plenum 

chamber 


8" Diam., 12' long 


| Pressure taps to manometers 




Air 

shutoff 

valve 


Sharp-edge 
orifice t 


8" Diam., 10.5' long 


4" Diam. 


Thermocouple 

rake Variable-area 
air orifice 
Dilution air 


4 d 

Laboratory 

air 

Exhaust 

control 

valves- 

4" diam. 


'DT 


Exhaust 

shutoff 

valve 


>4 shutoff 
valve 


T 

Laboratory 

exhaust 



Pressure- 
control 
butterfly 
valves- 
4" Diam 


Exhaust hood 
- -Free jet 
~ -Stub tunnel 
'-Spark electrodes 
--Turbulence-generating grid 

i r-J~~~20:l Contraction nozzle 

Fine -mesh screens 
Calming chamber 
Honeycomb 
Transition section 
Propane gas injection grid 
Thermocouple rake 


8" Diam., 15.5’ long 


K 


'■Explosion blowout diaphragm 
(b) 


Air- inlet 
shutoff 
valve 

<T/ 


(b) Free-jet installation. 

Figure I. — Air-ducting, wind-tunnel, and mean airflow 
instrumentation. 
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Air supplied from the laboratory air facility was metered 
with a standard variable-area orifice and controlled with 
the butterfly and gate valves shown in figure 1. Static 
pressures and temperatures were measured with conventional 
manometers and thermocouples. 

The ducting approach to the tunnel was designed from 
principles used in conventional low-turbulence wind tun- 
nels. The transition section was followed by a honeycomb 
and turbulence-damping screen in the calming section to 
reduce the turbulence present in the air-supply flow. The 
entrance to the enclosed tunnel (or free jet) was a nozzle 
with a contraction ratio of 20, and provisions were made 
at the tunnel throat for the insertion of turbulence-producing 
grids. Three grids were used for the investigation, all hav- 
ing a mesh to wire-diameter ratio of 5. The grid wire 
diameters were 0.0313, 0.063, and 0.125 inch. 

The enclosed-tunnel installation was equipped with an 
adjustable plug and nozzle by which critical flow could be 
maintained at the tunnel exit, so that the tunnel flow field 
could be isolated from pressure disturbances in the labora- 
tory exhaust facility. The free-jet installation was equipped 
with an exhaust hood into which diluting room air could be 
drawn in addition to the free-jet flow. 

fuel system 

The propane fuel used throughout the investigation had 
the following composition by liquid volume: 

Propane, percent 97. 76 

Ethane, percent - 

Isobutane, percent — 051 

The propane fuel system is shown in figure 2 and was iden- 
tical for both the enclosed-tunnel and free-jet installations. 
Propane was supplied from a laboratory facility and metered 
with standard rotameters. Pressure, temperature, and flow 
controls are shown in figure 2. The propane was injected 
into the airstream through a grid having 60 equally spaced, 
0 . 040 -inch-diameter holes. The injection grid was installed 
at a station upstream of the calming chamber as shown in 
figure 1. 

IGNITION SYSTEM 

The flame-globule ignition system consisted of a pair of 
spark electrodes protruding into the tunnel stream from 
opposing walls and an electrical energy source. The spark 
electrodes used in the enclosed tunnel were carefully stream- 
lined to a fineness ratio of 5 with a thickness tapering from 
0.019 inch at the tunnel center to 0.038 inch at the tunnel 
wall. The spark electrodes used in the free-jet installation 
were of circular cross section with a diameter tapering from 
0.015 inch at the stub tunnel center to 0.075 inch at the stub 
tunnel wall. The spark-gap spacing was 0.015 and 0.030 
inch for the enclosed-tunnel and free-jet installations, re- 
spectively. The spark-electrode positions for the two 
installations are shown in figure 1 . 

The spark-energy source was of the capacitive type and 
produced single sparks synchronized with the camera opera- 
tion in the enclosed-tunnel installation or successive sparks 
at the rate of 3 per second in the free-jet installation. The 
spark energy could be varied by changing capacitors. 


Propane injection gr id - 


r-|. 25 -In. pipe manifold 



Figure 2. — Fuel system. 


FLOW-FIELD INSTRUMENTATION 

In general, the flow-field instrumentation was adaptable 
to point measurements in either the enclosed-tunnel or free- 
jet streams. Mean-stream-velocity measurements were 
made with a conventional pitot-static probe and microma- 
nometer. Fuel-air ratio was measured directly by a mixture 
analyzer. Stream samples were drawn through a sharp-edged 
0 . 125 -inch-inside-diameter sampling probe to the analyzer, 
which indicated fuel-air ratio directly by the thermal- 
conductivity-bridge method. 

Turbulence measurements were made with constant- 
temperature hot-wire-anemometer equipment, which is 
described in appendix B. Sound pressure levels were meas- 
ured with a conventional microphone and a sound-pressure- 
level meter, also described in appendix B. 


free-flame-growth measurement 

Photographic instrumentation. — The growth of the free- 
flame globule was recorded photographically in the enclosed- 
tunnel installation with the optical system shown in figure 3. 
The photographic records were taken with a 16-millimeter 
Fastax camera operated at approximately 3600 frames per 
second and located at the end point of the parallel-beam 
schlieren system. The enclosed-tunnel windows were 
striation-free plate glass. The light source was a mercury- 
vapor BH-6 lamp operated with a high-voltage direct- 

current power supply. 
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Figure 3. — Photographic instrumentation for measurement of free- 
flame growth in enclosed-tunnel installation. 


Ionization-gap instrumentation. — The growth of the free- 
Hame globule was measured in the free-jet installation with 
an ionization-gap probe, a photomultiplier flame sensing 
unit, and two flame counters. The physical orientation of 
these instrumentation components is shown in figures 4(a) 
and (b). The photomultiplier unit and its accompanying 
flame counter recorded the total number of flames that 
occurred during any particular run. The ionization-gap 
probe and its accompanying flame counter recorded the 
number of flames that grew to a radius of r or greater at the 
axial station x during any particular run. 

The photomultiplier unit consisted of a 931 A photo- 
multiplier vacuum tube with a fixed-plate voltage supply. 
The ionization-gap probe is shown in detail in figures 4(c) 
and (d). The brass fairing along the probe tip prevented 




(a ) 


(a) Over-all view of installation. 
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f- 0 . 125 " 

*-0.0313" Inconel wire 


^-Magnesia in 
/ swaged probe 


kmaMBmmmiwiviivm im i jmitmi i. 


0.125- Outside -diam. 
Inconel tubing 

(d) 


- 0.250 -Outside-diam. 
Inconel tubing 


(c) Ionization-gap probe. 

(d) Construction details of ionization-gap probe. 


Figure 4. Ionization -gap instrumentation for measurement of free-flame growth in free-jet installation. 
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flame seating on the probe both by flame quenching and by 
elimination of the recirculation zone. Because of the con- 
voluted form of each flame globule, the voltage signal from 
the ionization gap usually contained multiple peaks, which 
caused spurious counts in the flame counter. In order to 
smooth these peaks, resistors, capacitors, and diodes were 
included in the 360-volt direct-current battery circuit for 
the ionization gap. 

In order to prevent spurious counts arising from small 
fluctuations in the voltage signal from both the photomulti- 
plier unit and the ionization-gap probe, a long time constant 
was included in the flame-counter design. These flame 
counters employed thyratron-controlled relays and solenoids 
which actuated a mechanical counter. 

Photomultiplier instrumentation. — Simultaneous measure- 
ments of the spark light intensity and ensuing free-flame- 
globule growth were made in the free-jet installation with 
the photomultiplier instrumentation shown in figure 5. The 
photomultiplier unit consisted of a 931A photomultiplier 
tube with an ad j ustable-plate voltage supply. Light emitted 
from the spark and ensuing flame passed through the narrow 
slit to the photomultiplier unit. The voltage signal from 
the photomultiplier unit was recorded on 35-millimeter film 
from an oscilloscope trace. The photomultiplier instrumen- 
tation and the assumption required to relate output voltage 
to flame speed are discussed in appendix C. 

PROCEDURE 

The test schedule was as follows : 


Grid, in. 

Stream 

velocity, 

u, 

ft/sec 

Propane-air 
weight ratio 

Distance 
downstream 
of grid, 
x, in. 

Wire 

diam. 

Mesh 

Free jet 

0. 125 
.125 
.063 
.0313 

0. 625 
.625 
. 313 
.156 

35 to 105 
70 
70 
70 

0.07 

. 056 to . 090 
. 058 to . 084 
. 053 to . 084 

15. 2 to 16. 2 
15. 2 to 16. 2 
15. 2 to 16. 2 
15. 2 to 16. 2 

Enclosed tunnel 

0. 125 

0.625 

35 to 140 

0.045 (approx.) 

6.3 to 15.3 


ENCLOSED-TUNNEL INSTALLATION 

Room air was supplied to the enclosed tunnel and ex- 
hausted at critical flow through an exit choke by means of 
the altitude exhaust facility. Tunnel static pressures and 
temperatures were measured by manometers and thermo- 
couple rakes as shown in figure 1 (a). Fuel was injected into 
the airstream and metered through a standard rotameter, 
after which dilution air was admitted downstream of the 
choke. The tunnel stream velocity was then measured witli 
a Pitot-static probe and micromanometer. 

The sequential schlieren photographs of the propagation 
of the free-flame globule in the fuel-air mixture were obtained 
by actuating a single switch which synchronized the following 
operations: fuel shutoff (producing a semi-infinite fuel slug), 


Spark 
electrodes 7 


r 1 Q03I3" Slit' 

j. 48"—— 



camera start, ignition by a single spark, and camera stop 
after an appropriate interval. 

free-jet installation 

An exhaust airflow rate of approximately 500 pounds per 
minute was initiated, and dilution air was admitted to the 
hood. Laboratory pressurized air and propane were ad- 
mitted into the stub tunnel, after which the fuel-air ratio 
was adjusted to the desired value by means of a probe and 
the'NACA mixture analyzer; the stream velocity was 
determined from Pitot-tube measurements. 

The ionization-gap probe was positioned at a station x 
downstream of the spark electrodes and at a radial distance r 
from the tunnel centerline. A series of flame globules was 
initiated by the ignition system, firing at the rate of 3 
sparks per second. The number of flame globules N P 
intercepting the ionization-gap probe and the total number of 
flame globules ignited N, were read from instrumentation 
shown in figure 4(b). 

FLOW-FIELD MEASUREMENTS 

In order to check the uniformity of the flow field, velocity 
profiles were measured in both the enclosed-tunnel and 
free-jet installations. The mean stream velocity for both 
installations was essentially constant over the core of 
interest, as shown by figure 6. Similar measurements were 
made to check the uniformity of the propane-air ratio in 
both installations. The propane-air ratio was found to be 
constant over the major portion of the enclosed-tunnel flow 
field, and typical profiles measured in the free-jet installation 
are shown in figure 7. 

In order to establish the characteristics of the turbulent 
field, extensive measurements of the turbulence intensity 
the spectrum function F, and the correlation coefficient 
j were made in both the enclosed-tunnel and the free-jet 
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Figure 7. Typical fuel-air-ratio profiles in free-jet installation. 




Figure 8. Longitudinal turbulence intensity measured in enclosed- 
tunnel installation. 

installations. The definitions of these quantities and the 
instrumentation used to measure them are discussed in 
appendix B. 

ENCLOSED-TUNNEL TURBULENCE MEASUREMENTS 

Longitudinal turbulence intensity measurements made in 
the enclosed-tunnel installation are shown in figure 8. The 
data shown in the figure agree well with the von Kdrm^n 
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(a) Enclosed-tunnel installation. Measured 10.75 inches downstream 
of grid of 0.125-inch-diameter wire by 0.625-inch mesh. 

(b) Free-jet installation. Reference velocity equals velocity at 
2/ = 0.5 inch and x=1.0 inch downstream of spark electrodes. 
Figure 6. — Typical mean-stream velocity profiles. 


PROPAGATION OF A FREE FLAME IN A TURBULENT GAS STREAM 


data quoted in reference 15 and show a substantial decrease 
in turbulence intensity from the electrode position (x/d=5 0) 
to the end of the flame-growth observation field (x/d= 122). 
Longitudinal intensity profiles measured in the tunnel 
centerplane perpendicular to the observation plane were 
essentially flat over the core of interest. Turbulent energy 
spectra were measured in the enclosed-tunnel installation as 
described in appendix B. The scale of turbulence was 
estimated by comparison of experimental data with a grid 
of spectrum density curves based on an exponential form of 
the longitudinal correlation coefficient /=exp (£/£/)• A 
typical measured spectrum is shown in figure 9, and the 
estimated turbulence scales are listed in table I. 

FREE-JET TURBULENCE MEASUREMENTS 

Conversion of the enclosed tunnel to a free jet resulted in 
an intense sound field originating in the dilution air and 
exhaust hood and piping. As described in reference 16, 
velocity fluctuations associated with sound waves contribute 
to the hot-wire-anemometer measurements through the 
relation . . 

with the assumption that no correlation exists between the 
turbulence and sound velocity fluctuations. In order to 
determine the applicability of equation (3) in the free jet, 
longitudinal velocity fluctuation intensities were measured 
both with and without the intense sound field. These meas- 
urements are shown in figure 10 for each of the three turbu- 
lence grids used in the investigation. The velocity fluctua- 
tion intensity was definitely increased with the presence of 
the sound field for each of the three turbulence-generating 


The sound pressure level in the free jet was measured with 
the instrumentation described in appendix B and was found 
to be 136 decibels, which corresponds to a sound velocity 



fluctuation V«f of 0.96 foot per second. Values of the sound 
velocity fluctuation were also calculated from the data 
of figure 10 by the use of equation (3). These values of Vwf 
are compared in table II with the value of V«! corresponding 
to the sound pressure level of 136 decibels; good agreement 
is shown in every case. 

The kinetic-energy spectrum of the sound field alone is 
shown in figure 1 1 ; it has a generally continuous form with 
peaks, or periodicities, at 57, 110, 140, 170, 230 to 285, 600 
to 870, 1050, 1500, 2400, 3000, 4000 to 5600, and 7000 to 
9000 cycles per second. Periodicities observed in the sound 
spectrum were generally found again in the longitudinal 
velocity fluctuation energy spectra measured in the free-jet 
installation, as shown by the typical spectrum in figure 12. 

Longitudinal double-velocity correlation coefficients ./ 
measured in the free-jet installation are shown in figure 13 
for each of the three turbulence-producing grids. The irregu- 
lar form of the correlation-coefficient curves is attributed to 
the periodic nature of a considerable portion of the velocity 
fluctuations. The correlation coefficients were measured by 
the special methods described in appendix B and are defined 
by the relation 

(4) 
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Figure 10. — Longitudinal velocity fluctuation intensity in free-jet 

installation. 
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Figure 10.— Concluded. Longitudinal velocity fluctuation intensit 
in free-jet installation. 



Figure 11 . Longitudinal kinetic-energy spectrum of free-jet sound 

field. 



J igure 12. Typical longitudinal velocity fluctuation energy spectrum 
measured in free-jet installation. Grid: 0.125-inch-diameter wire, 
0.625-inch mesh ; distance downstream of spark electrodes, 4 inches ; 
mean stream velocity, 68.5 feet per second ; static pressure, 29 inches 
of mercury absolute; static temperature, 545° R. 
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where the symbols are defined by the following vector 
diagram: 

u 0 


The longitudinal turbulence scale is usually defined by the 
relation 

L/am Io fd * (5) 

Because of the periodic nature of the combined turbulence 
and sound fields, the longitudinal correlation coefficient did 
not approach zero within the range of measurement. There- 
fore, equation (5) could not be used to define the turbulence 
scale. Since the growth of the flame globule was observed 
over times comparable with £ of 0.1 inch, only the first 
portions of the correlation curves shown in figure 13 were 
considered relevant to the experiment. As a part of the 
calculations involved in the Scurlock-Grover analysis, values 
of the lateral scale of turbulence L t were calculated as. 
described in appendix D. These calculated values of L g are 
listed in the table in the RESULTS section. 

DETERMINATION OF FLAME SPEED 
fundamental considekations 


The relation between flame speed and free-flame-globule 
expansion rate is derived in reference 13 for a constant- 
pressure laminar free flame propagating in a quiescent fuel- 
air mixture. This relation was obtained by equating the 
mass flow of the unburned gases entering a flame front of 
infinitesimal thickness to the mass flow of the burned gases 
leaving the flame front. In a physical sense, equation (2) 
accounts for both the flame-front motion due to flame propa- 
gation and the motion due to thermal expansion of the burned 
gases. The following modification of this equation to account 
for a finite flaine-tront thickness 8 is presented in reference 14: 




(6) 


For this investigation, the thickness of the flame front was 
considered negligible compared with the radius r so that 
equation (2) was used to analyze the data. The values of 
Ptt , computed from R a and T a , were obtained from measure- 
ments of tunnel approach-stream flow. The values of 
PF , computed from R F and TV for each fuel-air ratio, were 
taken from references 17 and 18, respectively. The value 
of prlpa varied from 0.12 to 0.14 in the present investigation. 

A number of assumptions are necessary to the definition 
of effective turbulent flame speed as given by equation (2). 
The definition requires that — 

(1) The free-flame globule grow uniformly in all radial 
directions 

(2) The combustion proceed in a constant-pressure field 




(b) Grid: 0.063-inch-diameter wire, 0.313-inch mesh. 



(c) Grid: 0.0313-inch-diameter wire, 0.153-inch mesh. 

Figure 13. — Longitudinal double- velocity correlation coefficients 
measured in free-jet installation. Mean stream velocity, 68.5 feet 
per second. 



Mean radius of flame globule, 7j in. 
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(3) The flame globule be swept downstream at the stream 
velocity 

(4) No circumferential flow exist about the flame globule 

(5) Variations in spark energy not affect the globule 
propagation rate 

The validity of the assumptions was established by special 
tests described later in this section. 

PHOTOGRAPHIC METHOD 

Photographic data obtained in the enclosed tunnel were 
analyzed by projecting the 16-millimeter motion-picture 
film record of each flame-globule expansion onto a film 
viewer. A typical photographic record is shown in figure 
14(a). The outline of each globule was traced in order 
to permit planimetering its area and calculating a mean 
radius. The mean radius was plotted against time for each 
flame globule of the sequence, as shown in figure 14(b). 

Start 



Finish 


(a) Typical photographic record. Film speed, approximately 3600 
frames per second. 
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Figure 14 —Typical turbulent-flame-speed data from en closed-tunnel installation. 
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The slope of a straight line faired through such a series of 
experimental points is dr/dt. The effective turbulent flame 
speed was then calculated by multiplying dr/dt by the expan- 
sion factor pf/puj as shown in equation (2). 

From sets of 30 or more flame speeds obtained in this 
manner for identical fuel-air flow conditions, the cumulative 
probability of occurrence of flames having these flame speeds 
was calculated. In order to illustrate the statistical nature of 
these data, the turbulent flame speeds of a typical set of data 
have been plotted against cumulative probability of occur- 
rence in figure 14(c). Each datum point represents the 
effective flame speed of a single free flame. Cumulative 
probability of occurrence is defined as the percentage of 
flames having flame speeds greater than a given value. 


each of the two sets of data as shown in figure 16. The 
difference A r was obtained from the curves for the cumu- 
lative probabilities of 2, 20, 50, 80, and 98 percent. Equa- 
tion (2) was then employed to determine the flame speeds 
corresponding to the five probability percentages. 

PHOTOMULTIPLIER METHOD 

The photomultiplier method was used to obtain values 
proportional to the effective flame speed. The method 
essentially consists of measuring the free-flame-globule 
growth rate along the axial diameter of the globule. The 
method is described in the APPARATUS section and in 
appendix C. The instrumentation is shown in figure 5. 

VERIFICATION OF EXPERIMENTAL METHOD 


IONIZATION-GAP METHOD 

In order to obtain a set of turbulent flame speeds for a 
particular condition of turbulent approach-stream intensity 
and fuel-air ratio in the free-jet installation, the ionization- 
gap probe was positioned consecutively at two stations 
z=15.2 and 16.2 inches, as outlined diagrammatically in 
figure 15. For each chosen position r, at least 1000 flames 
N t were ignited. The ionization-gap counter recorded the 
number of flames N P intercepting the probe. The ratio 
Np/N t was the cumulative probability of a flame globule 
having a radius at least as great as that indicated by the 
probe setting r. The ratio N P /N t was then plotted on prob- 
ability' coordinates against r and a curve was faired through 



The discussion of equation (2) in a preceding section in- 
cludes a number of assumptions necessary to the derivation 
of the equation. Since equation (2) essentially expresses 
the experimental method used in the investigation, the 
assumptions were verified by tests that are described in 
the following paragraphs. 

Uniformity of globule growth. — Photographic data of the 
flame-globule growth in the enclosed tunnel were analyzed 
to determine whether the globules had preferential growth 
directions. This analysis was made by measuring globule 



(a) Fuel-air ratio, 0.056; mean stream velocity, 68.3 feet per second. 


(b) Fuel-air ratio, 0.070; mean stream velocity, 69.6 feet per second. 


\CL ) X UUi“an v .ysw , i**vw** w w ~ A 

Figure 16.— Typical cumulative probability of flame-globule size measured with ionization-gap instrumentation in free-jet installation. Grid: 
0.125-inch-diameter wire, 0.625-inch mesh; static pressure, 29.3 inches of mercury absolute; static temperature, 545 R. 
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diameters (see fig. 17) at four angles of rotation from the 
tunnel axis. While each globule had different growth 
rates along the four diameters, no consistent relation be- 
tween growth rate and direction was found. This obser- 
vation supports the use of a mean radius in the photographic 
method of flame-speed measurement, the use of a single 
radial direction in the ionization-gap method, and the use 
of a single diameter in the photomultiplier method of flame- 
growth observation. 

Pressure pulsations. — In order to determine whether the 
fluctuations in flame-globule-growth rates, illustrated by 
figure 14(b), were due to pressure fluctuations in the en- 
closed tunnel, data as shown in figure 17 were inspected for 
possible phase relations between the fluctuations in growth 
rate along the four diametral directions. Since no frequency- 
phase relation was found from these data, it was concluded 
that pressure fluctuations were not responsible for the 
fluctuations in globule-growth rate. Additional measure- 
ments with a wall-mounted pressure pickup indicated the 
absence of significant pressure pulsations in the enclosed 
tunnel. 
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Figure 17. — Typical flame-globule-diameter histories measured in 
en closed-tunnel installation . 


Axial velocity of flame globule. — The effect of buoyancy 
of the hot flame globule on its axial velocity was determined 
by measuring the distance from the spark electrodes to the 
globule center on each frame of photographic sequences 
such as shown in figure 14(a). The axial velocities calcu- 
lated from such measurements agreed closely in all cases 
with the measured stream velocity, which indicated that 
the buoyancy effect was negligible. 

Unburned-gas motion. — As shown by equation (2), the 
rate of thermal expansion of the flame globule is from 7 to 
8 times as great as the effective turbulent flame speed. 
In order to determine whether the unburned-gas motion 
caused by this thermal expansion was purely radial (from 
the globule center), simultaneous oscillograph records were 
made in the free-jet installation with the ionization-gap 
probe and a shielded hot-wire-anemometer probe as shown 
by the typical oscillograph and probe position diagram in 
figure 18. The ionization-gap trace indicated the flame- 
globule passage, and the anemometer trace was propor- 
tional to the instantaneous stream velocity. Inspection of 
60 such records indicated no axial change in mean stream 
velocity and therefore no circumferential flow of unburned 
gases about the sphere. The marked decrease in amplitude 
of velocity fluctuations at the time of flame-globule passage 
is interpreted as being due to the radial translation of the 
free-jet stream core by the thermal expansion of the flame 
globule. 

Effect of spark energy. — In order to determine the effect 
of spark energy on the flame-globule growth, two series of 
runs were made in the free-jet installation under identical 
stream and fuel conditions but with spark source energies 
differing by a factor of 2.5. As shown in figure 19, cumu- 
lative probabilities of flame-globule sizes were measured at 
a series of axial positions with the ionization-gap instrumen- 
tation for the two spark source capacitors of 0.0004 and 
0.0010 microfarads. The data of figure 19 show an increase 
in initial flame-globule size, but no change in turbulent 
flame speed between the axial positions #=15.2 and 16.2 
inches. 

Effect of spark light intensity. — Additional information on 
the relation between the spark characteristics and the flame- 
globule growth was obtained with the photomultiplier in- 


lonization- 

gap 

i ' probe,-- 



Ionization-gap trace 


Time •* 

Figure 18. — Typical oscillograph of ionization-gap and hot-wire 
signals in free-jet stream at time of flame-globule passage. Grid: 
0.125-inch-diameter wire, 0.625-inch mesh; stream velocity, 70 feet 
per second; static pressure, 29 inches of mercury absolute; tempera- 
ture, 80° F. 
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strumentation shown in figure 5. With this instrumentation, 
oscillographs were obtained in a free-jet installation show- 
ing spark light intensity and flame light intensity as a 
function of time for 151 consecutive flame globules at con- 
stant mean stream and fuel conditions. Typical oscillo- 
graphs obtained by this method are shown in figure 20. 
As shown in appendix C, the portion of the trace due to 
the spark light intensity was proportional to spark current, 
and statistical analysis of the magnitudes of the peak spark 
light intensity indicated that this parameter approximated 
a normal probability distribution. The flame speed of the 
globule was taken as proportional to the slope of the portion 
of the oscillograph trace due to flame light intensity, as 
discussed in appendix C. 

The group of 151 consecutive runs is shown in figure 21, 
in which peak spark light intensity is plotted against flame 
speed for each particular flame globule. The statistical cor- 
relation coefficient for this data was determined with Pear- 
son’s product moment formula (ref. 19) and was found to 
have an absolute value of | 0.17 | for a possible range of zero 
to | 1.0 |. The low value of the statistical correlation coeffi- 
cient is interpreted to indicate that variations in spark 
current had no effect on turbulent flame speed. 

Additional verification of statistical variation of turbulent 
flame propagation. — The cumulative probability of occur- 
rence of slopes measured from the series of 151 photo- 



Figure 19. — Effect of spark source energy on flame-globule growth. 


multiplier oscillographs is shown in figure 22. The statis- 
tical distribution of flame speeds measured by this method 
is interpreted as additional proof of the randomness of the 
phenomena. 

Validity of common-origin assumption, — The ionization- 
gap method, described earlier in this section, assumes that 
all flame globules have a common origin point, as shown in 
figure 15. In order to determine the error introduced by 
this assumption, “origin points” were measured from the 
series of 151 photomultiplier traces by extrapolating the 
slope of the flame light intensity r trace to the time axis, as 
shown in figure 20. The origin points showed a statistical 
deviation in position, and the theoretical error introduced 
into flame speeds determined by the ionization-gap method 
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Figure 20. — Typical oscillographs of photomultiplier probe trace 
depicting spark occurrence and growth of flame globule in free-jet 
stream. 
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Figure 21. — Correlation between spark light intensity and turbulent 
flame speed for 151 consecutive flame globules at identical mean 
stream and fuel conditions. 
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was calculated to be ±2.5 percent or less, within one standard 
deviation of origin positions. The error in flame speed due 
to shift in origin point is small compared with the deviations 
in flame speed actually measured by the photographic 
(fig. 14(c)), the ionization-gap, and the photomultiplier-tube 
methods (fig. 22). 


PRECISION OF DATA 

Ionization-gap measurements. — Each datum point plotted 
as cumulative probability against flame radius, such as shown 
in figure 16, represents information from at least 1000 flames. 
As shown in reference 20, samples consisting of 1000 flames 
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Figure 22. — Statistical variation of flame speed as determined from 
photomultiplier oscillograms. 


would have the following accuracy in the cumulative prob- 
ability value: 


Cumulative 

probability, 

percent 

Possible 

error, 

percent 

50 

±5 

20 

±7 

2 

±29 


The possible errors listed correspond to 1.28 standard devia- 
tions of the error distribution and therefore should hold 80 
percent of the time. The general consistency of the data 
throughout the lower probability region such as shown in 
figure 16 seems to justify greater confidence in the low 
cumulative probability range than that shown in the preced- 
ing table. 

The reliability of the curves faired through the plotted 
experimental data, such as presented in figure 16, may be 
estimated from the deviation of data points from the fitted 
curves. 

The greatest divergence was found in results for the 0.125- 
inch-diameter wire. For the 98-percent probability level, 
the maximum data scatter from a faired curve was 20 per- 
cent; for the 50-percent level, 7.8 percent; and for the 2- 
percent level, 3 percent. The maximum scatter for the 
0.063-inch-diameter wire was 5.3 percent at 98-percent 
probability, and for the 0.0313-inch-diameter wire, it was 13 
percent at the 98-percent level. 

As discussed earlier in this section, the value of flame 
temperature used in the equation for turbulent flame speed 
was a calculated theoretical value. It may be of interest to 
note that flame temperature has a much greater effect on the 
laminar flame speed than on the calculation leading to 
effective turbulent flame speed. As shown by reference 
21, a decrease of theoretical adiabatic flame temperature from 
4000° to 3900° R would result in a laminar flame speed 
decrease of 43 percent. In contrast, this change in flame 
temperature would cause only a 2.5-percent decrease in the 
turbulent flame speed calculated from equation (2). It 
appears, then, for the purpose of this investigation, that the 
effect of a small error in actual flame temperature is neg- 
ligible in calculating the density ratio. 

Other sources of inaccuracy may be found in the instru- 
mentation and equipment. The positioning of the ionization- 
gap probe involved an error not exceeding 0.5 percent. 
Measurement of the fuel-air ratio was accurate within 2 
percent. Flow measurements were made with a precision 
of 3 percent. For the measurements characterizing the 
turbulence parameters, an estimated value of ±5 percent is 
suggested for the probable limit of error. 

Photomultiplier method. — A treatment of the sources of 
error present in flame-speed measurement by the photo- 
multiplier method is presented in appendix C. 

Photographic method. — In addition to the errors intro- 
duced by the equipment and instrumentation, sources of 
error peculiar to the photographic method must be con- 
sidered. Such factors as adjustment of the projected image 
to one-to-one correspondence, timing, optical resolution, and 


PROPAGATION OF A FREE FLAME IN A TURBULENT GAS STREAM 


15 


scale factor are all considered significant and have been 
thoroughly treated in reference 22, which assigns a relative 
maximum error of 10 percent to the photographic-reduction 
technique. 

RESULTS 

Effective turbulent flame speeds were measured in the 
enclosed-tunnel and free-jet installations over a range of 
propane-air and turbulence conditions as shown in the 
following table: 


Grid, in. 

Stream 
velocity, 
U, ft/sec 

Longitudinal 

velocity 

fluctuation 

intensity, 

Vw2, ft/sec 

Lateral 
turbulence 
scale, L t , in. 

Propane-air 
weight ratio 

Distance 
downstream 
of grid, x, in. 

Diam. 

Mesh 

Free jet 

0.125 

0. 625 

35 to 105 

1. 18 to 2. 60 

0. 0645 

0. 07 

15. 2 to 16. 2 

.125 

.625 

70 

1.69 

.0645 

. 056 to . 09 

15. 2 to 16. 2 

.063 

. 313 

70 

1.34 

.0587 

. 058 to . 084 

15. 2 to 16. 2 

.0313 

. 156 

70 

1. 17 

.0437 

. 053 to . 084 

15. 2 to 16. 2 

Enclosed tunnel 

0. 125 

0. 625 

35 to 140 

0. 81 to 3. 51 

0. 08 to 0. 13 

0. 045 

6. 3 to 15. 3 






(approx.) 



The lateral turbulence scales L g for the enclosed tunnel 
were obtained from the values of L f shown in table I by 
using equation (Blc) of appendix B. The values of L g for 
the free jet were obtained by fitting an exponential curve 
to the lateral correlation coefficient g as described in ap- 
pendix D. Both methods for evaluation of L g are based 
on the assumption of isotropy, which was not proven because 
of the lack of turbulence measurements in the lateral 
direction. 

The effective turbulent flame speeds, plotted in figure 
23, were measured in the free-jet installation at a fixed 
propane-air ratio of 0.07 over a range of mean stream 
velocities with the 0.125-inch wire diameter by 0.625-inch 
mesh turbulence- generating grid. The aata are plotted as 
a function of longitudinal velocity fluctuation intensity 
Vt? with cumulative probability as a parameter. Since 
turbulence measurements were made at only one stream 
velocity, and since the sound pressure level was inde- 


TABLE I— SUMMARY OF LONGITUDINAL TURBULENCE 

SCALES 

[Estimated from longitudinal energy spectra in enclosed-tunnel installation. Mean stream 
velocity, 50 ft/sec; static pressure, 26 in. Hg abs; static temperature, 540° R. Comparative 
values calculated through the isotropic turbulence relation Lr=2L t from lateral scale of turbu- 
lence data reported in ref. 15.1 


Distance down- 
stream of grid, 
x, in. 

Longitudinal scale of 
turbulence, Lf, in. 

Present in- 
vestigation 

Ref. 15 

6.3 

0. 16 


9.3 

.21 


12.3 

.25 

0. 17 

15.3 

.25 

.21 


pendent of stream velocity, the data of table II and the 
following equation were employed to calculate the longi- 
tudinal velocity fluctuation intensity: 

Vw*= -y/[y (^)] +«I (7) 

The turbulence intensity (^v%/U) was assumed constant 
throughout the range of stream velocity U (see ref. 15), 
and u* corresponded to the measured sound pressure level of 
136 decibels. 



Figure 23. — Turbulent flame speed as a function of longitudinal ve- 
locity fluctuation intensity with cumulative probability of occur- 
rence as a parameter. Propane-air weight ratio, 0.070; stream static 
temperature, 545° R; static pressure, 29.3 inches of mercury abso- 
lute; grid: 0. 125-inch-diameter wire, 0.625-inch mesh; free-jet 

installation. 


TABLE II.— COMPARISON OF VELOCITY FLUCTUATIONS IN TURBULENT STREAM 
WITH AND WITHOUT INTENSE SOUND FIELD 


[Velocity fluctuations measured along duct axis. Measured value of sound pressure level, 136 db; stream velocity, 68.5 ft/sec.] 


Grid, in. 

Distance 
downstream 
of turbulence - 
producing 
grid, x, in. 

Stream velocity fluctuations, 
ft/sec 

Calculated values of sound 
field 

Measured values of 
sound field 

Diam. 

Mesh 

With sound 
field, \ u 2 

Without sound 
field, 

V*df=V^-ttr 

1 'ft/sec 

Sound pres- 
sure level, 
db 

Ff 

VMi) 

ft/sec 

Sound 
pressure 
level, db 

0. 125 

0. 625 

3 

1. 81 

1.40 

1. 15 

137.5 

0.96 

136 



4 

1.75 

1.32 

1.15 

137.5 

.96 

136 

0.063 

0. 313 

3 

1.32 

0. 926 

0.94 

135.5 

0.96 

136 



4 

1.41 

.892 

1. 10 

137.0 

.96 

136 

0. 0313 

0. 156 

3 

1.12 

0.688 

0. 89 

135.0 

0.96 

136 



4 

1.22 

.619 

1.05 

136.5 

.96 

136 
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Figure 24. — Turbulent flame speed as a function of turbulence in- 
tensity. Propane-air weight ratio, approximately 0.045; stream 
temperature, 540° R; static pressure, 26 inches of mercury absolute; 
en closed-tunnel installation; grid: 0.125-inch-diameter wire, 0.625- 
inch mesh. 

Effective turbulent flame speeds measured in the enclosed- 
tunnel installation at a fixed propane-air ratio of approxi- 
mately 0.045 are plotted in figure 24 as a function of the 
longitudinal velocity fluctuation intensity Vw 2 with cumu- 
lative probability as a parameter. The data show a con- 
sistent increase of effective turbulent flame speed with 
velocity fluctuation intensity. These data were obtained 
downstream of the 0. 125-inch-diameter-wire by 0.625- inch- 
mesh grid over a range of stream velocities from 30 to 140 
feet per second. As shown by the data of figure 8, the 
decay of the turbulence intensity was quite severe over the 
range of interest (x/d=5 0 to 122). A turbulence intensity 
of (^Ut/U) =0.027 was taken as a representative value and 
assumed constant to obtain values of Vw 2 over the stream- 
velocity range. 

The effects of propane-air ratio and turbulence scale on 
the propagation rate of the free flame were investigated in 
the free-jet installation by measuring the turbulent flame 
speed downstream of the three turbulence-generating grids 
over a range of propane-air ratio from 0.053 to 0.090. (The 
variations in turbulence scale are listed in the table at the 
beginning of this section.) These data are shown in figure 
25 along with a dashed line representing laminar flame-speed 
data reported in reference 21. Turbulent flame speeds 
from the enclosed tunnel are included in figure 25(a) and 
were evaluated by interpolating the data of figure 24 at a 
velocity fluctuation intensity equal to that of the free-jet 
data in figure 25(a). 

(a) Grid: 0.125-inch-diameter wire, 0.625-inch mesh. 

(b) Grid: 0.063-inch-diameter wire, 0.313-inch mesh. 

(c) Grid: 0.0313-inch-diameter wire, 0.156-inch mesh. 

Figure 25. — Turbulent flame speed as function of fuel-air ratio with 

cumulative probability of occurrence as parameter. Mean stream 
velocity, 70 feet per second; stream static pressure, 1 atmosphere; 
static temperature, 85° F. 



.050 .055 .060 .065 .070 .075 .080 .085 


Fuel-air weight ratio 
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DISCUSSION 

MEASUREMENTS OF PRESENT INVESTIGATION 

A discussion of the effect of turbulence on flame propaga- 
tion is generally conceded to rest on the principles of turbu- 
lent motion set forth by Taylor in reference 23. These 
principles are summarized in the equation relating the mean 
displacement X of a fluid particle from its original position 
to two parameters of the turbulent field : 


where D is the turbulent diffusion coefficient in the in- 
direction, Xis the particle displacement in the indirect ion, and 
& is the Lagrangian double-velocity correlation coefficient 
defined by the relation 




u 0 u t 
m ' u 2 0 Via? 


( 9 ) 


It is of particular importance to note that the turbulence 
parameters in equation (8) are both mean values, and that 
instantaneous, local particle displacements depend on the 
statistical deviations about those mean values. It is sug- 
gested that the fluid particle displacements have a statistical 
distribution about the mean much the same as the statistical 
variation of velocity fluctuations in the turbulence. Such 
variations in particle displacements would be particularly 
noticeable when observed over times less than the char- 
acteristic time of the turbulence, as is the case in the present 
investigation. 

A statistical variation in the particle displacement X 
would probably cause variations in the speed of flames where 
flame wrinkling was controlling; a statistical distribution of 
the turbulent diffusion coefficient D would also probably 
cause variations in the flame speed where the diffusion process 
affects turbulent flame propagation. As a consequence of 



(a) Free-jet data. Propane-air ratio, 0.070. 


this statistical distribution of turbulence velocity and cor- 
relation coefficients, particular free flames might encounter 
conditions especially favorable and therefore have high 
propagation rates, while others might be subjected to condi- 
tions that would reduce their rate of growth. This picture 
of statistical variation is consistent with the data of the 
present investigation shown in figures 23 to 25 and also agrees 
with observations of stabilized flames that show a wide 
brush of flame fronts. 

The effect of flame temperature on laminar flame speed has 
been reported in references 21 and 24, wherein it is shown 
that small changes in flame temperature result in large 
changes in laminar flame speed. In particular, reference 24 
shows that this dependence of laminar flame speed on flame 
temperature is of greater magnitude in the lean and rich 
regions. The lower portions of the flame-speed probability 
band in figure 25 generally show decreases in the ratio of the 
turbulent to the laminar flame speed S T /S L in the lean and 
rich fuel-air regions. This decrease is consistent with the 
effect of flame temperature if the lower portion of the proba- 
bility band is considered to represent flames that have en- 
countered turbulent conditions sufficiently violent to de- 
crease the flame temperature by diffusive action. 

Comparison of figures 25(a), (b), and (c) shows that 
greater portions of the flame-speed band fall below the 
laminar flame speed as the grid size becomes smaller. As 



(b) Enclosed-tunnel data. Propane-air ratio, approximately 0.045. 


Figure 26. — Comparison of turbulent-flame-speed data with those of other investigations. 
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shown by figure 13, decreases in the grid size reduce the 
mean turbulent eddy size. Since the effect of the turbulence 
is thought to be increasingly diffusive in character as the tur- 
bulence scale decreases, this trend might be due to an in- 
creasing proportion of flames having temperatures less than 
the adiabatic value as the turbulence scale becomes smaller. 

COMPARISON WITH DATA FROM OTHER INVESTIGATIONS 

A comparison of the data of the present investigation with 
those of references 4, 8, 9, and 22 is presented in figure 26. 
In order to obtain a basis for comparison, the data are plotted 
with the ratio of turbulent to laminar flame speed S T /S L 
against longitudinal velocity fluctuation intensity V?. 
The turbulent flame speeds measured in a Mache (nozzle) 
type burner (ref. 4), a Bunsen burner (ref. 8), and a flame 
tube (ref. 9) fall above those measured by the free-flame 
method in references 14 and 22 and in the present investiga- 
tion. It should be noted that the S T /S L data from the 
present investigation are plotted against the longitudinal 
velocity fluctuation intensity, which includes the contribu- 
tions from both turbulence and sound. A better agreement 
with the data of other investigations could be obtained bv 
plotting the ratio S T /S L against the velocity fluctuation in- 
tensity Vu| due to turbulence alone. This procedure is 
questionable, since the sound velocity fluctuations may 
have an appreciable effect on turbulent flame speed when 
coupled with turbulence, and also because the magnitude of 
the sound field in the other investigations is unknown. 
The abscissas of figure 26 may be converted to turbulence 
intensity by the use of equation (3) and the data of table II. 

The effect of fuel-air ratio on turbulent flame speed of free 
flames is shown in figure 27, in which the data from the free- 
flame experiment reported in reference 22 are compared with 
those of the present investigation. The data from the present 
investigation shown in figure 27 were taken 15.2 to 16.2 
inches downstream of the 0.125-inch-diameter-wire by 0.625- 
inch-mesh grid, while the data from reference 22 were taken 
14 to 18 inches downstream of a 0.125-inch-diameter-wire by 
0.500-inch-mesh grid. The grid sizes and axial distances are 
close enough so that the turbulence scales were roughly equal. 
Because of a large difference in stream velocity, the two sets 
of data were taken under widely different velocity fluctuation 
intensities. Both sets of data show a general trend of de- 
creasing turbulent- to laminar-flame-speed ratio in the lean 
and rich fuel-air-ratio regions. The considerable deviation 
from mean values exhibited by the data of reference 22 sub- 
stantiates the statistical deviation of turbulent flame speed 
noted in the present investigation. 

COMPARISON OF DATA WITH VARIOUS THEORIES 

The current prevailing concept of turbulent flame propaga- 
tion is based on references 25 and 26, wherein the proposals 
are made that (1) small-scale turbulence (of the order of the 
flame-front thickness) increases the local flame propagation 
rate because of increased diffusion rates, and (2) large-scale 
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Figure 27. — Comparison of turbulent-flame-speed data with those 

of reference 22. 


turbulence increases the over-all flame propagation rate by 
wrinkling the flame with attendant increased flame-front area. 
Although it is well known that a turbulent field contains a 
statistical distribution of eddy sizes, the complexity of the 
problem has dictated the assumption of either the iarge- or 
small-scale case. 

The basic concept proposed by Shelkin (ref. 26) for the 
large-scale case is that the wrinkling of the flame front might 
be approximated by cones having a base proportional to the 
turbulence scale and a height proportional to the turbulence 
intensity. This concept has been supported experimentally 
by Hottel and co workers in reference 1. 

Turbulent flame speeds measured in flames stabilized on 
bluff bodies in high-velocity fuel-air streams led Scurlock and 
coworkers to postulate the concept of flame-generated turbu- 
lence (refs. 6 and 7). Theoretical analyses to account for 
flame-generated turbulence were proposed by Karlovitz and 
coworkers (ref. 2) and by Scurlock and Grover (ref. 11), but 
neither analysis conforms to all existing experimental data. 
A recent analysis by Tucker (ref. 27) predicts that little, if 
any, turbulence would be generated by free flames with small 
surface deformations. 

Since the Shelkin equation is not only of historical interest 
but is also used as a basis for the Scurlock-Grover theory, 
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the data of the present investigation are compared with the 
Shelkin equation in figure 28, wherein the ratio S T /S L is 
plotted as a function of longitudinal velocity fluctuation 
intensity for four cumulative probabilities. The data include 
flame speeds measured with each of the three turbulence- 
generating grids and are all for the constant fuel-air weight 
ratio of 0.07 (figs. 23 and 25). The Shelkin equation was cor- 
related with the data at each value of cumulative probability 
by the method of least squares. This correlation was accom- 
plished by determining the constant B in the Shelkin equation 



The value of the constant B was estimated to be unity by 
Shelkin but varies from 0.867 for the cumulative probability 
of 2 percent to 0.220 for the cumulative probability of 80 
percent. The value of the constant B is shown in figure 28 



(a) Cumulative probability of occurrence, 2 percent. 

(b) Cumulative probability of occurrence, 20 percent. 


for each of the cumulative probabilities. The data scatter 
about the Shelkin equation but show some agreement with 
the Shelkin theory. A curve plotted through the data ob- 
tained with the 0.125-inch-diametcr-wire grid does not follow 
the shape of the Shelkin curve, however. 

The theoretical expression for turbulent fiaine speed pre- 
sented by Tucker in reference 27 was obtained by a wave 
analysis of the effect of a low-intensity turbulent field on an 
infinitesimally thin flame front having convolutions, or wrin- 
kles, of low slope. The final equation shows no effect of 
turbulence correlation coefficient and is, in the present 
notation, 

%-l+<Kr)« O') 

Oz, O/, 

where C{t) is a function of the ratio r= T F /T a . The analysis 
includes the assumption that the second-order terms ( u 2 /S‘l ) 2 
are negligible, that the flame temperature and laminar flame 
speed are dependent only on fuel-air ratio, and that the diffu- 
sive action of turbulence is negligible. 

In order to obtain a comparison with the present data, the 
Tucker analysis was extrapolated well beyond its intended 
limit of the ratio (u?/S L ). The data of the present investiga- 
tion compare well with the Tucker analysis in the region 
of low velocity fluctuation intensity, as shown in figure 29, 
but fall below the theory in the higher-intensity range, where 
the theory cannot be expected to hold. A further comparison 
of the turbulent-flame-speed data with the Tucker analysis 
is made in figure 30, wherein the ratio S T /S L is plotted 
against fuel-air weight ratio. These data are for each of the 



(c) Cumulative probability of occurrence, 50 percent. 

(d) Cumulative probability of occurrence, 80 percent. 


Figure 28. — Comparison of turbulent-flame-speed data with theory of reference 26. Fuel-air weight ratio, 0.07. 
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(b) Grid: 0.063-inch-diameter wire, 0.313-inch mesh. 


(a) Grid: 0.1 25-inch-diameter wire, 0.625-inch mesh. 


Figure 30.— Comparison of turbulent-flame-speed data with theories from references 11 and 27. 
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Figure 29. — Comparison of turbulent-flame-speed data with theories 
of references 11 and 27. Propane-air weight ratio, 0.07; stream 
static pressure, 29.3 inches of mercury absolute; temperature, 
545° R; grid: 0.125-inch-diameter wire, 0.625-inch mesh. 


three turbulence-generating grids and were taken at the same 
mean stream velocity. The data in each figure fall below 
the theory, but, with the exception of the enclosed-tunnel 
data (fig. 30(a), fuel-air ratio, 0.045), follow the general 
trend of the Tucker analysis in the upper portion of the 
cumulative probability band. These data suggest that the 
analysis might be taken as the limiting value of turbulent 


flame speed for the turbulence conditions specified in the 
figures. 

The turbulent-flame-speed data are also compared with 
the Scurlock-Grover analysis (ref. 11) in figures 29 and 30. 
This analysis assumes that the flames proceed locally at the 
laminar rate S L} but that the over-all flame speed S T is 
increased because of flame-front wrinkling caused by tur- 
bulence. The geometrical model is similar to Shel kin’s, 
except that the wrinkle height is assumed to be positively 
affected by turbulent motion and negatively affected by 
local flame propagation. The analysis also accounts for 
flame-generated turbulence caused by velocity and density 
gradients across the flame front. 

Flame-generated turbulence has been neglected in a 
modified Scurlock-Grover analysis, which was used in the 
present investigation and is discussed in appendix D. The 
measured correlation coefficients, as shown in figure 13, were 
used in the analysis, in contrast to the exponential forms 
assumed in the Scurlock-Grover report (ref. 11) and the 
comparison with free-flamc data made in reference 22. The 
analysis predicted that the turbulent flame speed had not 
yet reached its full value at the axial stations of measure- 
ment (x=15.2 and 16.2 in.). Therefore, values of S T /S L 
were calculated from the theory for flame-growth times 
corresponding to both stations and are shown in figures 29 
and 30. 
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(c) Grid: 0.0313-inch-diameter wire, 0.156-inch mesh. 

Figure 30. — Concluded. Comparison of turbulent-flame-speed data 
with theories from references 1 1 and 27. 


The data of the present investigation fall below the values 
of S t /S l predicted by the modified Scurlock-Grover analysis. 
The data follow the Scurlock-Grover analysis more closely 
for the largest turbulence scale (fig. 30(a)), and the agreement 


becomes progressively worse as the turbulence scale becomes 
smaller (figs. 30(b) and (c)). This type of relation may be 
due to an increasing diffusive action of the turbulence as the 
scale becomes smaller, as discussed previously in this section. 

The data agree somewhat better with the Scurlock-Grover 
analysis than with the Tucker theory in the extremes of the 
lean and rich fuel-air-ratio regions. This probably is due to 
the increasing V u 2 /S L in those regions that tend to exceed 
the limits of the assumptions made in the Tucker analysis 
discussed in a preceding paragraph. 

In summary, values from both the Tucker and the 
Scurlock-Grover analysis seem to form an upper limit to the 
data of the present investigation, except in the extremes of 
the lean and rich fuel-air-ratio regions. 

CONCLUSIONS 

From the results of the investigation, the following conclu- 
sions are made : 

1. The propagation rate of a free flame in a turbulent field 
with intense random sound velocity fluctuations is a statis- 
tical variable having a probability of occurrence band 
bounded by a zero propagation rate and possibly an upper 
propagation rate well above the laminar flame speed. 

2. The propagation rate of a free turbulent flame is 
substantially less than those previously measured for sta- 
bilized flames. 

3. The statistical spread of free-flame propagation rates 
is markedly wider in both the lean and rich fuel-air-ratio 
regions. 

4. Values calculated from the Tucker analysis for free- 
flame propagation rates apparently form an upper limit for 
the statistical spread of free-flame propagation rates in the 
regions where the velocity fluctuation intensity is low. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, May 7, 1956 
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APPENDIX A 

SYMBOLS 


A cross-sectional area 

A smoothed time-average surface area of turbulent 

(lame 

7> constant in Shclkin equation 

(J constant 

D turbulent diffusion coefficient, sq ft/sec 

d wire diameter of turbulence-producing grid, in. 

E (1-c voltage, v 

VF rms of a-c complex wave voltage, v 
F longitudinal energy-spectrum density function, sec 

/ longitudinal double-velocity correlation coefficient 

<j lateral double-velocity correlation coefficient 

i hot-wire current at stream conditions, amp 

K gain of decade amplifier 

L f longitudinal scale of turbulence, in. 

L g lateral scale of turbulence, in. 

Lagrangian scale of turbulence, in. 

N number of flame globules 

n frequency, cps 

R gas constant 

Lagrangian double-velocity correlation coefficient 
r radial distance from tunnel axis, in. 

S flame speed, ft/sec 

SPL sound pressure level, db 

T temperature, °R 

t time, sec 

U mean stream velocity, ft/sec 

u longitudinal velocity fluctuation, ft/sec 

v,w lateral velocity fluctuations, ft/sec 

X fluid particle displacement due to turbulent motion, 

ft 


x distance downstream of turbulence-producing grid, 

in. 

y distance from tunnel wall, in. 

a angle, deg 

13 power band width of wave analyzer, cps 

5 flame-front thickness, ft 

£ spatial separation of velocity vectors comprising 

correlation coefficient, in. 

P density, lb/cu ft 

a one standard deviation 

r ratio of flame temperature to unburned-gas tempera- 

ture, TpfTa 

<p equivalence ratio, fuel-air ratio/stoichiometric fuel-air 

ratio 

12 hot-wire resistance at stream conditions, ohms 

Subscripts: 

a ambient, or approach stream 

F flame 

/ pertaining to correlation coefficient/ 

(j pertaining to correlation coefficient g 

L laminar 

n pertaining to a frequency interval 

o no-flow condition in anemometry, or reference point 

in space or time double-velocity correlation 
coefficients 

P pertaining to number of flames of a given radius or 
greater 
s sound 

T turbulent 

t total, or time 

£ at a distance £ from zero position 


APPENDIX B 


TURBULENCE AND SOUND INSTRUMENTATION AND ANALYSIS 


THEORETICAL BACKGROUND 


Turbulent fields are generally defined by the two param- 
eters intensity and correlation coefficient. The turbulence 
intensity is defined as the rms value of the turbulent velocity 
fluctuations VP, and V W 2 - The double-velocity corre- 
lation coefficient may be represented by a second-order 
tensor which for a homogeneous, isotropic field of turbulence 
can be reduced to two components as shown in reference 28. 
These two components are the longitudinal and lateral 
double-velocity correlation coefficients defined as follows: 


Longitudinal: f— u °^ 

u 2 



(4) 


Lateral: g=^r 

v 2 

where / and g are the longitudinal and the lateral correlation 
coefficients, respectively, u is turbulent velocity fluctuation, 
£ is the spatial separation between the points of measure- 
ment of the u, and the bar signifies a long time average. 

Longitudinal and lateral turbulence scales have been 
defined by the relations 

;,=/> ( 5 ) 

r* oo 

4=J o S<% 



(5a) 
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These scales are intended to represent a mean turbulent 
eddy size, and for homogeneous, isotropic turbulence may 
be related through the equation (ref. 28) 


,+i 

8 


(Bla) 


(Bib) 


which, by integration, becomes 

For the commonly observed exponential form of /=exp 
equation (Bib) becomes 

L f =2L g (Blc) 

Taylor has shown (ref. 29) that the one-dimensional turbu- 
lence energy spectral density function F is the Fourier 
transform of the longitudinal correlation coefficient/ 




w 2irnx , 
t cos —jj - dn 


(B2) 


where F is defined by the relations 


and 


/» 00 

5 -J„ 

j: 


F dn 
F dn= 1 


where u* represents the total kinetic energy in the longitudinal 
direction, and v\ represents the kinetic energy contained in 
the frequency interval n to n+dn. 

HOT-WIRE ANEMOMETRY 

The constant-temperature hot-wire anemometry equip- 
ment is shown diagrammatically in figure 31, and for the 
most part has been described in references 16 and 30. The 
equation relating the rms voltage E \ measured by the 
average square computer to the turbulence intensity has 
been discussed in reference 16 and is as follows: 



4 iU 

(2K)n(i 2 -i 2 0 ) 




(B3) 


The spectrum density function F was measured from the 
tape-recorder playback signal by the following equation 
which is given in reference 16: 



(B4) 


where E* is the square of the rms voltage reading of the 
wave analyzer. 

The double-velocity correlation coefficient/was determined 
from a double-channel magnetic-tape recording of the 


Soundproofing 



Figure 31. — Hot-wire anemometer and auxiliary instrumentation. 


instantaneous voltage signal from the hot-wire anemometer. 
The double-channel tape was played back into a double 
correlation computer which electronically performed the 
instantaneous multiplication of velocity fluctuations indi- 
cated in equation (4). The variable separation £ between 
velocity vectors was attained by the use of a movable and a 
fixed playback head. The separation £ between the velocity 
vectors was equal to the separation between the heads 
multiplied by the ratio of the stream velocity to the tape 
speed. The equivalence of the correlation coefficient 
measured by this method and by the two-probe method was 
proven as a part of a research program reported in reference 
31. 


SOUND-FIELD INSTRUMENTATION 

The sound instrumentation consisted of a standard 
microphone and sound-level meter from which sound-pressure 
level could be determined. The microphone signal was also 
impressed on the wave analyzer to obtain the sound kinetic 
energy spectrum through equation (B4). The relation 
between sound pressure level SPL and sound velocity 
fluctuations is given in reference 32 (for simple harmonic 
waves) as 


SPL=20 log (paV u 0 + 7 4 (Bo) 
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APPENDIX C 

PHOTOMULTIPLIER INSTRUMENTATION 


RELATION BETWEEN FLAME LIGHT INTENSITY AND FLAME SPEED 

Assumptions. — The conditions assumed in relating flame 
light intensity as indicated by the voltage output of a 
photomultiplier tube to the surface area of a flame are 
treated in reference 33. They may be summarized as 
follows: 

(1) It is postulated that, if a point source of light is in- 
creased to a finite size but the distance from the light source 
to the detector is constant and large in comparison with the 
maximum dimension of the flame, the light flux impinging 
on the detector will obey the inverse square law. Under this 
condition, a change in the surface area of the flame will pro- 
duce an approximately linear change in the flux intensity 
registered by the detector. 

(2) The flame must consist of a homogeneous cloud of 
nonabsorbing emitters. 

If the propagation of a frcc-flame globule moving in a gas 
stream at stream velocity is observed by a detector through 
a long slit placed axially with respect to globule travel, some 
additional requirements are: 

(3) The detector must have an essential^ flat light response 
over the entire field of view. 

(4) The slit must be narrow enough in relation to flame- 
globule diameter to approximate a diameter of the projected 
spheroid. 

(5) The flame-globule center must not suffer any appreci- 
able radial displacement from the tunnel centerline. 

(6) The flame globule must not have a preferential growth 
direction. 

(7) The flame temperature must be constant for all the 
flame globules. 

Conditions (1) and (4) were met by suitable placement and 
dimensioning of the apparatus. For condition (3), it was 
determined that the photomultiplier unit had a maximum 
light response error of 5 percent over the total light accept- 
ance angle. Photographic evidence obtained in the pre- 
liminary phase of the investigation indicated satisfactory 
fulfillment of requirements (5) and (6). Assumption (7) 
has been previously treated in the section Measurements 
of Present Investigation. 

Description of apparatus. — A 93lA-type photomultiplier 
having an adjustable-plate voltage power supply was 
mounted in a light-tight box behind a 0.25-inch-diameter 


aperture. This detector unit was positioned in an enclosed 
tunnel 4 feet from the free-jet stub tunnel, as shown in 
figure 5. The detector aperture was centrally located with 
respect to the slit in the machined aluminum plate, which 
replaced an upper portion of one wall of the stub tunnel. 
The voltage output of the photomultiplier was amplified 
and fed to an oscilloscope. Oscillograms were taken by 
high-speed, continuous motion-picture technique. 

Oscillograms. — As shown in figure 20, the oscillograph 
trace of the ignition and growth of a single flame consists of 
two portions, a spark light amplitude decay curve and a 
flame light amplitude growth curve whose slope at any point 
is equal to the change of diameter of the flame globule with 
time, governed by the considerations listed previously. 
This slope 2 dr/dt is proportional to turbulent flame speed 
S T as given by equation (2). 

RELATION BETWEEN SPARK LIGHT INTENSITY AND SPARK CURRENT 

Standring and Looms (ref. 34), using short-duration 
sparks produced by a lightly damped oscillatory current, 
have shown that the light emitted by a given spark gap 
depends on the current in the gap and its time variation. 

In order to determine the effect of a long-duration, direct- 
current pulsed spark on the current-spark light relation, an 
auxiliary investigation was undertaken in collaboration with 
Clyde C. Swett, Jr. The spark production and recording 
apparatus is fully described in reference 35. The photo- 
multiplier detector apparatus was as described previously 
except that the detector unit was mounted at one end of a 
rectangular tube 2 feet in length. The field of view of the 
photomultiplier tube was limited by a 0.0625- by 2-inch 
aperture in the other end of the tube. A Polaroid Land 
Camera was used to make oscillograms. 

A series of oscillograms was taken of sparks ignited in still 
air and in an airstream having a velocity of 50 feet per seeond. 
Spark duration was obout 600 microseconds, spark energy 
was of the order of 17 milli joules, and peak spark current 
approximated 0.23 amperes. A comparison of spark light 
amplitude-time curves with the current-time curves obtained 
in the apparatus designed by Swett showed a marked simi- 
larity of shape between the sets. It was concluded from this 
information that, at least for a first-order approximation, 
spark light intensity was roughly proportional to spark 
current. 
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APPENDIX D 


CALCULATION OF FLAME SPEED FROM SCURLOCK-GROVER ANALYSIS 


The basic equation in the Scurlock-Grover theory (ref. 11) 
is obtained by assuming that the increase of flame speed is 
due to the wrinkling of the flame front in the form of right 
circular cones, so that 



(I>1) 


where -^X 2 lL\ is assumed proportional to the ratio of cone 
height to base diameter and C x = 4 is a constant. The mean 
square particle displacement is assumed to depend on three 
effects : 

(1) Turbulent mixing, which increases the flame-front area 

(2) Laminar flame propagation, which decreases the flame- 
front. area 

(3) Flame-generated turbulence, which increases the tur- 
bulent mixing 

The first effect on X 2 is given by 

(D2) 

at Jo 


where Si is the Lagrangian correlation coefficient and g t is 
the Eulerian correlation coefficient. The second effect is 
given by 


)]"■-.} <w» 


d-Jx l_ t 

dt. 


where C*=Z\ is a constant. The combination of the first 
two effects is then 


d^X 2 1 dXj <hjxi 

dt -y/m dt dt 


(D4) 


Since the Lagrangian correlation coefficient cannot be meas- 
ured directly, Scurlock and Grover assumed it to be of the 
form given in reference 1 5 : 

^=exp^~ 


where jSf is the Lagrangian scale of turbulence and was as- 
sumed to be equal to % L e . In the present analysis, was 
assumed to be related to / through the independent variable 


transformation |=0.6 -yju 2 t (ref. 36). Isotropy was as- 
sumed, so that g was obtained from / through the relation 
(ref. 28) 


q=f+- — 
0 


(D6) 


The correlation coefficient g t was then obtained from the inde- 
pendent variable transformation %=S L t. The turbulence 
scale L g was calculated by matching the curve exp (— Z/L g ), 
g=e-i /L ‘, In g=L g i to the g plot previously calculated 
from /by equation (D6). 

All integrations and differentiations were done graphically, 

and the calculation procedure followed the following order: 

(1) Calculate g from equation (D6). 

(2) Transform/® to &[£/( O.Cy-Ju 2 )] 

(3) Transform </(£) to g,(£/S L ). 

(4) Calculate 


dt 


(5) Calculate 


2 u 2 j St g t dt 


(6) Calculate 




(7) Calculate 

d^T 2 _ 1 dX T d-yjTj 

dt dt dt 


(8) Calculate 

(9) Calculate 




dt 



As mentioned in the DISCUSSION section, flame-gener- 
ated turbulence has been assumed negligible, and the final 
form of the Scurlock-Grover theoretical equation is as given 
previously. 
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